GravitoMagnetic Force in Modified Newtonian Dynamics 

Qasem Exirifard 
School of Physics, Institute for Research m Fundamental Sciences, Tehran 19538-33511, /rariJ 

We derive the gravitomagnetic field in the ACDM and Modified Newtonian Dynamics (MOND) 
paradigms before showing that the the gravitomagnetic force at the edge of a galaxy can be in accord 
with only one of them. We notice that MOND enhances the precession of a gyroscope embedded in 
the gravitational saddle point of the Sun- Jupiter system, and the near future technology can detect 
the enhancement. 

PACS numbers: 04.20.-q,04.50.Kd,95.35.+d 
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The missing mass problem in galaxies can be resolved 
either by modifying the Newtonian Dynamics and grav- 
ity, or assuming the presence of dark matter. Perhaps 
it is interesting to search for some effects/predictions of 
MOND which can not be reproduced by the dark matter 
paradigm. In this line of thought, we show that MOND 
not only affects the fall of of the gravitational potential 
but it also affects the fall of of the gravitomagnetic force. 
We report that the ACDM and MOND paradigms pre- 
dict different fall of for the gravitomagnetic force on and 
beyond the edge of a spinning galaxy. We notice that 
MOND also alters the gravitomagnetic force near the 
gravitational saddle 'points' of the solar system: at the 
events where the gravitational field vanishes. We argue 
that the MOND paradigm can be approved or discarded 
in the near future by measuring the gravitomagnetic field 
in the saddle 'point' of the Sun- Jupiter system. 

The note is organized as follows. First we review the 
gravitomagnetic force before uplifting it into the MOND 
paradigm. We then estimate the size of the gravitomag- 
netic force in the saddle point of the Sun- Jupiter system. 
This suggests that the gravitomagnetic force in this semi- 
Lagrange point can be measured in the near future. This 
suggestion adds to the possibility of observing MOND in 
the solar system [l|, y] ■ 



I. GRAVITOMAGNETIC FIELD IN MOND 

Classical gravity is governed by a single scalar field, 
the gravitational potential. The Newtonian gravitational 
potential satisfies: 



V^* = 47rGp, 



(1) 



where p is the density of matter. Albert Einstein at- 
tempting to uplift gravity to a relativistic regime, first 
replaced the space-time metric of Minkowski by 



ds^ = ~c{<^fdt^ + dx^ + dy^ + dz^ 



(2) 



later with the Gromann's help, he introduced the Rie- 
mannian metric. 



ds' 



g^.^dx^^dx", 



(3) 



as the relativistic gravity 3]. The relativistic theory of 
gravity has a symmetric rank-two tensor: the metric. 
Metric has 10 components in four dimensions, 9 more 
than the degrees of the classical gravity. To perceive 
the physical meaning of the degrees of the freedom of 
the relativistic theory, let the trajectory of a slow mov- 
ing particle be considered in a static deviation from the 
Minkiowki metric. In so doing, the metric reads 



ds" 



ri^i,dx^dx^ -f h^j_vdx^dx^ , 



(4) 



Also note that due to the equations of motion of the 
metric, the off-diagonal components of hij are compara- 
ble to its other components only for a rela,tivistic mass 
distribution, mass distributions like geons ^]. Also note 
that the contribution of the hij can not be neglected on 
light propagation or orbits of massless particles. Here 
we, however, are considering the geometry around a non- 
relativistic mass distribution. We also study the orbits 
of massive particles: the paths that at the leading and 
sub-leading order can be derived from the first variation 
of 



(5) 



S = / rfr((-l + 2(t))t' + {iy + 2Aiti') 



wherein hijX^x^ has been ignored, and r is an affine pa- 
rameter and 



2^100 



(6) 

A., = hM (7) 

The Euler-Lagrange equation for t derived from ([5]) reads 

t(-l + 2Aq) + 2^ii' = cte (8) 

where i{—\ + 2^40) stands for the gravitational redshift 
while 2AiX^ represents a new relativistic term, a term yet 
to be experimentally observed [2j|. 

The Euler-Lagrange equation for x^ derived from ^ 
then leads to 



2^i 
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d'x 

dt^ 



c'd'Ao + cS''{dkA^~d,Ak) 



dx 

~dr 



(9) 



Now let it be defined 



The effective action then reads 



2$, 

500 = -(l + ^-j 



9oi 



2A, 
c 



(10) 

(11) 



using which the equation ^ can be rewritten as follows 



-V$ + w X (V X A) 



(12) 



This allows interpreting V x A as a gravitomagnetic field. 
Now consider a gyroscope in the space time. V x A cause 
precessions of the orbits of a test particle. This precession 
is referred to as as the Lense-Thirring precession Q ■ Ref. 
[8| provides a decent recent review on Lense-Thirring pre- 
cession for planets, satellites in the Solar system. 

Now consider a gyroscope in an orbit. The similar- 
ity between the gravitomagnetic field and magnetic field 
beside the spin precession formula in electrodynamics 
(S = fi X B,fi — -^S) dictates that its spin precesses 
byi 



n 



LT 



--V y.A 
2 



(13) 



This precession is called the Pugh-Schiff frame-dragging 
precession [iflliil- The Pugh-Schiff frame-dragging pre- 
cession due to the rotation of the earth recently has been 
measured by the gravity probe B with the precision of 
19% [H. 

In the linearized Einstein-Hilbert gravity, the field 
equations in the harmonic gauge simplifies to 



d^U'" = 



-IGTrGT^o; 



where hij is the trace reversed perturbation 
The linearized equations can be derived from 



(14) 
(15) 



(16) 



S 



d\{-]^dji^,.d''h^"' + le^G^r^" + Kd^Ji''") 



(17) 
where \i, is a local Lagrange multiplier enforcing (llSp . 
and T^'' is understood to represent the linear energy- 
momentum tensor. Note that the effective action is in- 
variant under the residual symmetry of the harmonic 
gauge. It is invariant under 



h^,u -^ h^^ + 5^^!. + 9,.Cm 


(18a) 


□?p = 


(18b) 



We do not fix the residual symmetry. We consider it as 
the symmetry of the action. We decompose hfj,^ as what 
follows 

h^vdx'^dx'' = -Aodt^ + Aidx'dt + hijdx^'dx' (19) 



j d^x{Cg+Cs+Cx) 



C, = 327rGA^ J'' + 167rG/iyT*^' 



where 






d^A, 



d.A^ 



(20a) 

(20b) 

(20c) 
(20d) 



(21) 
(22) 



Note that the effective action is invariant under the resid- 
ual symmetry of the harmonic gauge pS|) . The t/(l) 
part of the action - the part that includes A^ - resembles 
the ordinary electrodynamics. It is called the Gravito- 
ElectroMagnetism (GEM). Compared to the electrody- 
namics, the equations of motion for A impose two extra 
conditions on A,, 



9^Ap - 
doA,-djh{^ = 



(23) 
(24) 



The first one states that the GEM should be solved in 
the Lorentz gauge. The second condition ((24)) implies 
that GEM has wave solutions only if hij field possesses 
a wave solution. The wave solution of the theory is due 
to the dynamics of hij field. This means that though 
the GEM is akin to the ordinary electrodynamics, due 
to its constraints and in the absence of the tensorial per- 
turbation, it does not exhibit the wave solutions of the 
ordinary electrodynamics. 

Near and around the galaxies, hij is suppressed due to 
the non-relativistic velocity of the stars and gas inside the 
galaxy. At the leading order hij does not affect the orbits 
of slow-moving massive particles. Slow moving particles 
see only the GEM part of the metric (O. Since we are 
interested in the orbits of slow moving massive particles 
around a galaxy we thus consider only the GEM part of 
mi: 



d'^xiCA + Xd^A^) 



Ca = -^F^.F^" 



327rGA^J^ 



(25) 
(26) 



A time dependent A^ , through the constraint equation 
(|24p , induces a time-dependent behavior for hij . The or- 
bits of the stars are blind to the change in hij . In the 
study of the orbits of the stars, therefore, the time de- 
pendent solutions of ([25]) are valid. The reader should 
be cautious that gravity waves can not be addressed in 
the truncated action: waves require considering the full 
action (pO)) . The truncated action ([25l) is, however, suffi- 
cient for studying the orbits of the stars in galaxies. The 
symmetry of the truncated Lagrangian (j26p is 



A., 



A., 



d^A 



(27) 



where A is a general scalar field. 

The Newtonian approximation of the linearized GEM 
action (1^ reads 



A^ix,t) = (4$(a:),0) 



(28) 
(29) 



The effective action in the Newtonian approximation 
then simplifies to 



5 = 



d^a 



1 



SttG" 



|V$| 



p$ 



(30) 



In the Modified gravity reafization of the MOND [13J, 
one replaces the Newtonian classical field theory with a 
general field theory but retain the Newtonian dynamics 
{F = ma): 



SmoG — — 



d^a 



1 



SvrG 



/:($,v$,---) + P* 



, (31) 



Keeping intact the Newtonian dynamics means that or- 
bits of slow moving particles are derived from ([5]) . The 
AQUAL approach :14] assumes that the symmetries for 
the equations of motions derived from S and SmoG are 
the same. The symmetries for S are 



$ ^ $ + A. 



(32) 



where A is constant. Imposing (|32|) on ((3T|) requires C 
to be a functional of the derivative of the Newtonian 
potential: 



C = £(V$,V2$, 



(33) 



AQUAL also requires the equations to be second order. 
So the Lagrangian is simplified to 



C = /:(v$) 



(34) 



We can construct only one scalar out of V<I>. So the 
Lagrangian reads 



C = J"(V$.V$) 
and the AQUAL action follows 

1 



S AQUAL = 



d^a 



SvrG 



J-(V$.V$/a^) + p$ 



(35) 



(36) 



The first variation of the AQUAL action with respect to 
$ yields 



where 



V"(^(l ^) V„$) =47tGp, 



fj,{x) — J-'{x ). 



(37) 



(38) 



The MOND terminology than requires [16j [2; 



H{x) : 
p{x) 



1 : For X > 1 
! X : For x < 1 



and 



ao = (1.0 ±0.2) X 10^ 



^10' 



(39) 
(40) 

(41) 



In the following we would like to apply the AQUAL 
assumptions, as reviewed above, on ([20| . In this Letter 
we are, however, interested in the orbits of slow moving 
particles. The slow-moving particles feels A^ component 
of the metric ([T9| . What slow- moving particles test can 
be derived from (P5|) . The symmetries of ([25)) is (|27l) . 
In the line of AQUAL model, we search for a non-linear 
generalization of (|25|) leading to second order differen- 
tial equations. This generalization must coincide to the 
AQUAL model for vanishing gravitomagnetic field. 

The simplest non-linear Lagrangian density for A pre- 
serving (|27l) follows 



-MOND 



c 






327rGA^J'' (42) 



which must converts to (jSHj) in the absence of the gravito- 
magnetic field. The consistency with the AQUAL model 
(|5H|, therefore, requires 



C 



MOND 



= J" 






327rGA^ J^ (43) 



Near a static solution representing a slow rotating mass 
distribution where |VAo| ^ |V x A^j, we can taylor ex- 
pand (02]) around the gravitoelectric field strength and 
keep the leading and the subleasing terms: 

/:a/Ond(V$,VxA) « P(V$.V$/a2) 

IV$I 
+/i(^^)|VxAp (44) 
ao 

in which psp is utilized. Now let the equations of motion 
for ((31]) be computed 



-Vx 



|V$| 

ao 

V<I)| 

Ijl{- ^)V X A 

ao 



47rGp 
IGttG- 



-J 



(45) 
(46) 



where in (|45l) . terms proportional to the gravitomag- 
netic field strength squared have been ignored. In the 
'Newtonian' regime (15^ . the Einstein-Hilbcrt GEM are 
recovered. In the MOND regime both of the electric 
and magnetic field strengths are altered. We can write 
the electric and magnetic gravitational field strengths in 
the MOND theory in terms of their counterparts in the 
Einsten-Hilbert gravity: 



^^{ 



lv$| 



ao 



)V$ = V<^EH + V X ft, 



m( 



ao 



)v X A = V X Aeh + yh 



(47) 
(48) 



where h and h are identified by 



= V X V$ 



V.V X A 



^ V£^H + Vx^ (49) 



We could obtain (|47)) and (|48)) using an effective field 
theory approach. Let us illustrate this. Since the orbits 
of slow moving particles are sensitive to the gravitoelec- 
tric and gravitomagnetic fields, one can assume that their 
effective Lagrangian reads 



^MGA 



(Pr 



inG 



/:($,V$,--- ,A,dA,---) 



(51) 



which in the Newtonian regime must coincide to 



5' = - 



d^r 



1 



8ttG 



|V$p + |V X A\^)+p<^ + Aj.A 



(52) 

where is the gravitoelectric potential and A represents 
the gravitomagnetic potential. Assuming that ([51]) and 
([52]) have the same symmetries, one concludes that (|5T|) 
possesses the following symmetries 



$ -^ $ + V X A. 
A ^ A + VA, 
(Aj) ^ (-!,-/) 



These symmetries require 

£($,V$,--- ,A,dA,- 



^{Eg,Bg) 



(53) 
(54) 
(55) 



(56) 



where E = Vcj) and Bg = \7 xA. Note that we are assum- 
ing the same field content in the MOND and Newtonian 
regime. Theories likes TeVeS 17[ assume different field 
contents but they claim to have provided an interpolating 
theory from the very weak regime of gravity to the very 
strong regime of gravity. Here we adhere to the minimal 
assumptions: the only degrees of freedom in the weak 
regimes of gravity that can be observed are the pertur- 
bation around the metric. We, however, provide neither 
an interpolating theory to the strong gravity regimes nor 
the full generally covariant theory. 

Requiring ([5^)1 Lagrangian to coincides to ([55)1 leads 
to 

CiEg,Bg) « Ti\E^) + C^i^^)\Bg\^ 
Oq ao 

wherein the gravomagnetic field strength is assumed to 
be small: 



\Bg\ « \Eg\ 



(57) 



The equations of motion derived from (|57l) can be ap- 
proximated by 






V,E 



t^EH 



(58a) 



ao 



V X ( Ci{ \ ' )Bg 1 ^ VxBeh (58b) 

Note that due to ^} , we have used /^(■^^^-^) ~ K^)- 
We can use ([55)1 to express the fields in terms of the 
gravitoelectric and gravitomagnetic fields of the Einstein- 
Hilbert theory: 



^J■ 



Tp2 I n2 

-^ ^E = Eeh + V xh (59a) 



ao 



El + Bl 
Ci{^ '-)B = BKH + Vh 



ao 



(59b) 



where h and h are identified by the consistency conditions 
similar to (gH) and (|50)) . 

Now consider the y frame that is moving with the small 
velocity V to the x frame: y* = a;* -I- V^t. In this frame 

Eg = Eg + V xBg (60a) 

Bg = Bg-V X Eg (60b) 

Eeh = Eeh + VxBeh (60c) 

Beh = Beh-VxEeh (60d) 

We assume the same physics holds in all the inertial 
frames. So we have the same equations of motion in 
the X and y frame: 

Tp2 _|_ r>2 

m(^ -)Eg ^ Eeh + V x ^ (61a) 

ao 



ao 



-)Bg ^ BsH + Vh 



(61b) 



Inserting (|60| in (|6T|) . and using ([59l) leads to 



{ti-Ci)VxBg = Vx{h-h) (62a) 

{liL~Ci)VxEg = V{h-h) (62b) 

Recalling that Vi is a constant arbitrary small vector, the 
consistency for ([62]) demands 

(63) 
(64) 



SI.{^i-Cl)Eg = 
V X {fi-Ci)Bg = 



We demand the consistency for any given solutions, for 
arbitrary Eg and Bg. This requires 



^l^ Ci 



(65) 



Inserting which in (|59|) reproduces (|47|) and (j48|) . This 
perhaps independently affirms the correctness of our pro- 
cedure in obtaining (|T7)) and pS)) . 



II. GRAVITOMAGNETISM OF A SPHERICAL 
MASS DISTRIBUTION 

In the following we would like to compute the gravit- 
omagnetic and gravitoelectric fields that a slow rotating 
spherical mass distribution produces in the vacuum. In 
the Einstein-Hilbert gravity we have 



V<5eh 
V xAeh 



Gm _ 






(66) 
(67) 



where m is the total mass and J is the total angular ve- 
locity of the spherical mass, and f — represents the 
center of the mass distribution. Due to spherical symme- 
try V X h = satisfies (|49l) leading to 



|V$ 



Gm_ 



(68) 



which in MOND regime (|40|) can be explicitly solved for 



V* = (Gmao)2 — 



(69) 



This is the ordinary MOND modification of the Newto- 
nian field capable of resolving the missing mass problem 
in galaxies and reproducing the TuUy- Fisher relation 18]. 
In order to compute the gravitomagnetic field strength 
in the MOND regime around a slow rotating spherical 
mass distribution, we first compute V/i. Inserting (|69|) 
into the consistency equation for V/i, ([SO)). in the MOND 
regime (|40)) . then yields 



V.{TVh) = - — 



(70) 



which is the non-homogeneous Laplace's equation in four 
dimensions written in the spherical coordinates: 



nJ.^S^ 



(71) 



where h = h{r,r.J) is understood. Let it be emphasized 
that ([TT]) represents the h equation in large r. In the 
'Newtonian' regime, at the vicinity of the origin, it holds 
Vh = 0. So we choose the solution of dTlT) which is 
source free at the origin. Doing so, the fall of of the Vh 
is guaranteed to be r"^ or less. The gravitomagnetic field 
in the MOND regime, therefore, yields 



V X A 



G 



2c2(Gmao)2 r2 



3-/-5) + 0(^) (72) 



Note that ([72]) is not divergent for small masses because 
\J\ ex TO. It is ([7^ that describes the gravitomagnetic 
field strength around a spherical galaxy at its MOND 
regime. We note that the fall of of the gravitomagnetic 



field strengths of MOND is r~^ while that of the Einstein- 
Hilbert theory is r~^ . The gravitomagnetic field is en- 
hanced in the deep MOND regime (|40|) . 

The equations for the Newtonian potential and the 
gravitomagnetic field of the ACDM theory read 



V^* 



47rG(/9 + pDark) 






(73) 
(74) 



where /ODark and jDark are respectively the density and 
the angular velocity distributions of dark matter. The 
gravitomagnetic field strength that the ACDM theory 
predicts for a spherical spinning galaxy at its edge then 
follows 



V X A 



ACDM 



G_ J_ 

2c2S-3 



3^-^) 



(75) 



There exists no direct information available about the 
velocity distribution of dark matter. The theoretical sce- 
nario considers dark matter halo as a cloud under a ra- 
dial collapse. Such a collapse would not contribute to 
the angular momentum. We, additionally, observe that 
the difference between MOND and ACDM can not be 
simply assigned to the total angular momentum of the 
dark matter. We, therefore, conclude that measuring the 
gravitomagnetic force at the edge/beyond the edge of a 
galaxy refutes one of the MOND and dark paradigms 
and proves the other one. However the gravitomagnetic 
force at the edge of a galaxy is too small that one may 
not hope for its detection in the near future. In the next 
section, we argue that the gravitomagnetic force can be 
measured in the MONDian regime of the Solar system 
with the help of the near future technology. 



III. GRAVITOMAGNETIC FIELD IN THE 
SOLAR SYSTEM 

In some points in the solar system the gravitational 
field of the planets and the Sun and the galaxy cancel 
each other. Let these points be called the semi-Lagrange 
points. At a neighborhood around the semi-Lagrange 
points MOND/MOG governs gravity. Let these neigh- 
borhoods be called the MOND windows. If we send/keep 
a satellite to/in MOND windows we then can manifestly 
search for the prediction of the MOND paradigm by per- 
forming on board experiments. 

Perhaps one can identify the MOND windows of the 
Solar system and their properties by numerically solving 
(|47l) for the Solar system. We, however, need only to 
estimate the size of the MOND windows and the gravit- 
omagnetic field within them. The largest MOND window 
occurs at the semi-Lagrange point of the Sun- Jupiter sys- 
tem. In order to estimate the size of this window let us 
.just consider the gravitational field of Sun and Jupiter 
[231 . Assume that the Jupiter's orbit is a circle around 
the sun. The magnitude of the Newtonian gravitational 



field of tlie Sun and Jupiter in tlie line connecting these 
two reads: 



?(^) 



GMi 



Sun 



GMjupitor 



{d-xf 



(76) 



where x is the distance from the Jupiter and d is the 
distance between Sun and Jupiter. The Newtonian grav- 
itational acceleration vanishes at 



a;sL = 



1 



A/su 



(77) 



A/,, 



xsL identifies the Semi-Lagrange point. Around this 
point the Newtonian gravitational acceleration can be 
approximated to 



g{xsL + 5x)^-£l^^^^5x 



(78) 



or equivalently 



, , r N 2GMjupiter f Msu 
lixsL + Sx) « ( — 



d3 



'M 



Jupiter 



-f/^ Sx (79) 



where . f^^"° 

AjJupin 



::^ 1 is assumed. This states that the 
MOND window in the semi-Lagrange point of the Sun 
and Jupiter system is about 5km. This is large enough 
to accommodate a satellite. 

In the following we would like to estimate if we can 
measure the gravitomagnetic force produced at the Sun- 
Jupiter MOND window. We calculate this in the rest 
frame of the MOND window. The gravitomagnetic 
field receives contribution from the angular spins of the 
Jupiter and Sun, and also the orbital motion of the Sun 
and Jupiter. In the rest frame of the MOND window the 
leading contributions to the gravitomagnetric field are 
due to the orbital motion of the Sun and Jupiter : 



A 



Sun 



Aju 



pitcr 



^ G Afsun 2n 

^ G A/,Kipitcr 27T 



(80) 
(81) 



where T is the orbital period of the Jupiter around the 
Sun: 



where L represents the spin of the earth 

L w 7.1 X lO^^kg m^/s (84) 

and Tpf, is the radius of the orbit of the gravity probe B: 



rpb 



(6367.5 + 642)fcm ^ 7009.5km 



(85) 



So the leading contributions to the magnitude of the 
gravitomagnetic force at this MOND window yield 



A 



Sun 



^Probo B 
^Jupiter 
obo B 



Ap, 



2.1 X 10" 



7X 10" 



(86a) 
(86b) 



Note that these gravitomagnetic fields are due to the 
Lorentz boosts. Their contribution to the gyro scope pre- 
cession sometimes is attributed to Thomas [13| • They are 
relativistic effects. However since MOND enhances the 
gravitoelectric field in the MOND regimes, so it changes 
the Thomas precession. 

The magnitude of the gravitomagnetic fields ((86)) 
shows that utilizing the today technology we can not 
detect the gravitomagnetic field in the MOND window, 
even if we manage to place a satellite in its deep MON- 
Dian regime and enlarge the magnitude of the gravit- 
omagnetic field by some factors. The factor of 10^^, 
however, is not very tiny. Jupiter is not far from the 
earth too: there was a proposal ^] to measure the the 
Lense-Thirring effect caused by the spin of Jupiter in 
the ongoing JUNO mission f2J] , and there exists the new 
E JSM-Laplace mission to Jupiter [g^ . We might improve 
the current precision by a factor of 1000 in ten years pro- 
vided that we annually improve the precision by a factor 
of two. So in the near future perhaps we have reached the 
desired precision. When this precision is achieved then 
placing a gyroscope inside the MONDian regime of the 
gravitational saddle point of the Sun- Jupiter and measur- 
ing the Pugh-Schiff effect will prove or refute the MOND 
paradigm. 



T = 4331.6 days 



(82) 



Let us compare these contributions with the gravitomag- 
netic field measured by the gravity probe B. The gravity 
probe B measures the gravitomagnetic field around the 
earth through the Lense-Thirring effect. It has recently 
measured the gravitomagnetic field around the earth: 



G L 



(83) 



Acknowrledgements 

I thank ICTP for the nice hospitality during which this 
work was accomplished. I thank J. Bekenstein and L. 
lorio for the comments and criticisms through the corre- 
spondences , S. Randjbar-Daemi, M. Sheikh- Jabbari and 
H. Firouzjahi for the discussions. 



[1] Jacob Bekenstein and Joo Magueijo , Phys. Rev. D 73 
(2006) 103513 ; M. Milgrom, arXiv:0906.4817 



[astro-ph.CO]. 



[2] L. lorio, Open Astron. J. 3 (2010) 156 [arXiv:0901.3011 

[gr-qc]]- 

[3] H. F. M. Goenner, arXiv:0811.4529 [gr-qc]. 

[4] T. Jacobson, Class. Quant. Grav. 24 (2007) 5717 

[arXiv:0707.3222 [gr-qc]]. 
[5] John Archibald Wheeler, Phys. Rev. 97 (1955) 511. 
[6] R. Penrose, Rev. Nuovo Cimento, 1 (1969) (Special 

Number), 252. 
[7] J. Lense and H. Thirring, Phys. Zeits. 19, 156 (1918). 
[8] L. lorio, H. I. M. Lichtenegger, M. L. Ruggiero and 
C. Corda, Astrophys. Space Sci. 331 (2011) 351 
[arXiv: 1009.3225 [gr-qc]]. 
[9] T. Padmanabhan, Gravitation: Foundations and 
Frontiers, Cambridge University Press, 
www.cambridge.org/9780521882231 . 
[10] Pugh, G.E.: WSEG Research Memorandum No. 11 

(1959). 
[11] L. I. Schiff, Phys. Rev. Lett. 4, 215 (1960). 
[12] C. W. F. Everitt et al, Phys. Rev. Lett. 106 (2011) 

221101, arXiv: 1105.3456 [gr-qc]. 
[13] M. Milgrom, Astrophys. J. 270 (1983) 365; M. 

Milgrom, Astrophys. J. 270 (1983) 371. 
[14] J. Bekenstein, M. Milgrom, Astrophysical J. 286, 
(1984) 7. 



[15] B. Famaey and J. Binney, Mon. Not. Roy. Astron. Soc. 

363 (2005) 603 [arXiv:astro-ph/0506723]: 
[16] J. D. Bekenstein, arXiv:astro-ph/0701848; K. G. 

Begeman, A. H. Broeils and R. H. Sanders, Mon. Not. 

Roy. Astron. Soc. 249 (1991) 523. 
[17] J. D. Bekenstein, Phys. Rev. D 70 (2004) 083509 

[Erratum-ibid. D 71 (2005) 069901] 

[arXiv:astro-ph/0403694] . 
[18] R. B. Tully and J. R. Fisher, Astron. Astrophys. 54 

(1977) 661. 
[19] SL. H. Thomas, Phil. Mag. 3, 1 (1927). 
[20] Lorenzo lorio. New Astronomy 15 (2010) 554. 
[21] This is the kind of the modification of the effective 

energy of a particle that lets the extraction of energy 

from a black hole ( the Penrose mechanism) [y|. 
[22] The most widely used from /i are [11, [ill : fJ.{x) = ^^ 

and fi{x) = ^ 



[23] Including the gravitational fields of other planets would 
alter the position of the window but it is very unlikely 
to significantly change its size. 

[24] http://www.nasa.gov/mission_pages/juno 

[25] http://sci.esa.int/ejsm 



